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Abstract. A detailed treatment of the entanglement dynamics of two distant but 
non-identical systems is presented. We study the entanglement evolution of two 
remote atoms interacting independently with a cavity field, as in the double Jaynes- 
Cummings (J-C) model. The four-qubit pairwise concurrences are studied, allowing for 
asymmetric atom-cavity couplings and off-resonant ineractions. Counter to intuition, 
imperfect matching can prove advantageous to entanglement creation and evolution. 
For two types of initial entanglement, corresponding to spin correlated and anti- 
correlatcd Bell states $ and a full, periodic and directed transfer of entanglement 
into a specific qubit pair is possible, for resonant interactions, depending on the choice 
of relative couplings. Furthermore, entanglement transfer and sudden death (ESD) can 
be prevented using off-resonant interactions, although for some initial states, detunings 
will trigger an otherwise frozen entanglement, to allow a full entanglement transfer. 
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1. Introduction 

Entanglement is now regarded to be a resource central to the development of quantum 
technologies. It has been extensively studied theoretically and creation of entanglement 
has now been reported in a range of systems, including trapped ions p], atomic 
ensembles [2] and photon pairs [3]. 

A fundamental feature of entanglement is that it is easily degraded when the 
entangled system interacts with another [H EJ Ej. Work by Yu and Eberly [HI [7] has 
shown that the entanglement of two qubits can deteriorate rapidly, to the point of an 
abrupt total destruction, when coupled to an environment that results in irreversible loss. 
Recent experiments [H [9] confirm the sudden death of entanglement (ESD) and there 
have been further studies investigating this behaviour [TO], [TTJ [121 ESI EES HU El [16l [18] . 

Particularly relevant for the purpose of quantum networks are entangled qubits that 
can be stored, localised and controlled. Many investigations have focused on the storing 
of entanglement in two level atoms, that are constrained to intertact with a single mode 
field in a cavity [19]. The use of very high Q cavities and Rydberg atoms [20] minimises 
the irreversible losses of the type studied by Yu and Eberly. In this case, the atom 
interacts only with the single cavity mode and the Jaynes-Cummings model [2TJ, [22] is 
realized. The dynamics and control of entanglement between two isolated atoms A and B 
in a "lattice" , allowed only to interact with two "control" fields a and b, repsectively, as 
in the JC model is thus of fundamental interest [231 12H [25] . 

The evolution of the entanglement between two isolated atoms that are each coupled 
to a field via Jaynes-Cummings cavity has been studied recently by Yonac et al [231 121] 
and Sainz and Bjork [TJ]. By examining the pairwise entanglement using the method 
of concurrences [26], Yonac et al [23] showed that the atom-atom entanglement is 
transferred to the cavity-cavity system, and then back again, in a periodic fashion. Two 
types of initial two-qubit entanglement of the atoms were considered, the Bell state \l/ 
where the spins are anti-correlated, and the Bell state $ where spins are correlated. A 
complete and abrupt loss of atom-atom entanglement (ESD) is found to occur for $, 
provided the initial entanglement is not maximum, though the full entanglement is 
periodically regained. 

In this paper, we study the effect of non-identical JC systems, so that one can 
control the entanglement between A and B, via manipulation of the localised atom- field 
interaction or detunings. We present complete entanglement solutions, for the case of 
asymmetric atom-cavity couplings and off-resonant interactions (detunings), so that all 
possible pairwise concurrences are calculated. The results reveal that a full periodic 
and directed transfer of the entanglement from one qubit pair to another occurs with 
suitable choice of coupling ratios. We find that the full entanglement transfer and ESD 
can be controlled with the use of suitable detunings, for both initial states \1/ and <3>. 
Thus the general effect of detuning is to stabilize entanglement. For some initial states 
however, where entanglement is already constant, the effect of detuning can be to trigger 
a maximal entanglement. 
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2. Hamiltonian and energy states 

We begin by introducing the Hamiltonian of the system and the energy basis states that 
we will use to solve for the entanglement creation and evolution. 

2.1. The Double Jaynes-Cummings Hamiltonian 

We consider a system composed of two separated single mode cavities each containing 
a single two-level atom, as shown in Fig. [TJ In general, the two subsystems may not be 
identical in that the cavity frequencies and the coupling strengths of the atoms to the 
cavity modes could be different. 



Subsystem 1 Subsystem 2 




Cavity a Cavity b 



Figure 1. Schematic diagram of the system considered, two distant and non- 
interacting subsystems 1 and 2 each with a single-mode cavity containing a two-level 
atom. 

The Hamiltonian for the system in the electric-dipole and rotating-wave 
approximations can be written in the Jaynes Cummings (JC) form 

H = H Field + H Atom + H inU (1) 

where 

H Fie i d = hw x (a^a + ^)+hu; 2 (b ] b + (2) 

is the Hamiltonian of the cavity modes with resonant frequencies io\ and uj 2 respectively, 
H Atom = hio S A + h\JoSg (3) 

is the Hamiltonian of the atoms of the transition frequency ujq, and 

H int = hg x (atSj + aSi) + % 2 (b^S B + bSty (4) 

is the interaction Hamiltonian between the atoms and the cavity modes. The parameters 
gi and g 2 are the strengths of the coupling between the atoms and the cavity modes, that 
can be assumed to be real with no loss of generality. The and S~ are respectively 
the raising and lowering operators of the ith atom, and (a) and (6) are the creation 
(annihilation) operators for the modes of the cavities labelled a and b, respectively. 
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In previous treatments [231 [21], the focus of study has been on identical and resonant 
cavities, where U\ = u>2 = Wq and g\ = g2- This makes the two subsystems 1 and 2 non- 
distinguishable, which might suggest a greater robustness for entanglement creation. By 
contrast, we will show that asymmetric cavities and cavity-atom detunings can prove 
advantageous in enabling a control of entanglement. Our objective then is to include 
7^ ^2 7^ &a an d g\ 7^ g2, which also has the purpose of better modelling a real 
experimental situation, where it may be difficult to produce identical cavities. The 
unequal coupling constants for example may arise when atoms are not in equivalent 
positions inside the cavities. In addition, we find the case of asymmetric couplings leads 
to a breakdown of a simple entanglement transfer conservation rule derived by Yonac 
et al [21] , and proved to hold where g\ = . 

2.2. Energy basis states 

We assume an initial entanglement between the distant atomic systems A and B, and 
study the evolution of the entanglement in the absence of external losses, where there 
is a local coupling of each atom to its cavity mode. The energy basis states of the non- 
interacting Hamiltonian are the product states of those for the atoms and the cavity 
modes. We denote the energy states of a two-level atom as | j) and | j), for the excited 
and ground states respectively, so that S z \ f) = || f) and S z \ j) = — 1| J.). The energy 
states of a cavity are denoted by the Fock states |n), with n being the number of quanta 
in the cavity mode. The size of the Hilbert space for the system depends on the number 
of energy quanta involved in the evolution of the atoms and the cavity modes. 

We follow the approach of Yonac et al. [23], and consider two different types of 
initial entanglement, corresponding to the Bell state with spins anti-correlated 

\m) = cosa|ti 00) +e i/3 sina||t 00), (5) 

and the Bell state with spins correlated 

|$) = cosaltt 00) + e i/3 sina||| 00). (6) 

Here, | fl ij) represents the atom A in the excited ("up") state, the atom B in the 
ground ("down") state, and % and j photons in the cavity modes a and b, respectively. 
We will use the notation * and |\l/) interchangeably, to represent the quantum Bell 
state. Following Yonac et al. [23], we allow for an initial non-maximal entanglement 
state by incorporating the parameter a. 

Firstly, in Section [31 we consider the initial state |\1/) where there is only a single 
excitation present, so that the space of the system is spanned by four state vectors, 
defined as follows: 

1^) = |U 00), |H/ 2 ) = |IT00), |* 3 ) = mio), |tt 4 ) = |U01). (7) 

In addition, we will study the "delocalised" entangled states of type J2t=i l^)/2j where 
the component states have equal weightings, to show that such states can give a stable, 
albeit reduced, entanglement. 
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Secondly, in Section HJ we will consider the case where there are two excitations 
present, one in each system, so the space of the system is spanned by the vectors 

= |TT oo), |$ 2 ) = |noi), |$3> = UTio>, |$ 4 ) = |Un>. (8) 

In this case we also include the ground state of the system 

|0> = |$ > = I II 00) (9) 

for which there is no excitation present, so that we can specifically study the case where 
the initial entanglement is in the form of the Bell state (jSJ). 

We will quantify entanglement by using the measure of concurrence, introduced by 
Wootters [26]. The concurrence is defined as 

C = max (o, - \A2 - \As - V^l) , (10) 

where A« are the eigenvalues of the density matrix p' = p(a y ® o- y )p*(cr y ® a y ), and a y is 
the Pauli matrix in y direction. The maximum possible entanglement is given by C — 1, 
while C = implies separability. The state is entangled (inseparable) for C > 0. One 
can calculate that the concurrence for the states |\&) and |$) is 

C\i = 2 sin a cos a, (11) 

to indicate maximal entanglement for a = 7r/4. 

The transfer of entanglement between the two qubit subsystems is governed by 
rules derived recently by Yonac et al. [21], and Chan et al. [2T] . Sainz and Bjork [H] 
have also presented an entanglement invariant, for the case of the double JC model. 
Chan et al. [22] have proved the following result for entanglement initially in the form 
of the Bell state 

Cab + @Ab + + Cab = Ci2- (12) 

The total non-local pairwise entanglement, as measured by the concurrence squared, is 
conserved, where Cu is defined as in Eq. ( TIT]) to be the bipartite entanglement between 
the systems we label 1 and 2. In fact, from the work of Chan et al., it is seen that this 
rule applies to any initial state for which there is a single-excitation (refer to Section [3]). 
For the case of symmetric JC couplings, where g\ = gi, Yonac et al [21] have proved 
that in fact the following simpler result also holds, if the atoms are initially in the Bell 
state |^): 

Cab + C ab = C\i. (13) 

For the case of entanglement initially in the form of the Bell state |$), the pairwise 
entanglement is not conserved [2H], and following inequality is always valid [22] for 
non- interacting systems 1 and 2: 

< C\ B + C\ h + Cl B + Cl b < C? a . (14) 
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3. Evolution of single-excitation states 

3.1. Dynamical solutions 

We consider first the dynamics of the system in a pure state where only one excitation 
is present, so that the state of the system is always of the form 

\m T {t)) = cosa|10) + e i/3 sina|01). (15) 

The entanglement properties of such a state for the case of perfect coupling have been 
studied in detail by Yonac et al. [23l [24] . We allow for the possibility of imperfect 
couplings of the atoms to the cavity modes, where detunings can be non-zero and 
coupling constants different. For the purpose of solution, we introduce the wave function 
of the system 

|W(t)) = d 1 (t)\ n oo) + d 2 (t)\ it 00) + d 3 (t)\ I| 10) + di(t)\ II 01), (16) 

which is a linear combination of the available product states (0) of the atoms and the 
cavity modes at time t. The coefficient di(t) determines the probability amplitude of 
the ith state at time t. 

We find the coefficients di(t) by solving the Schrdinger equation 

l hj t m)) = H\m- (17) 

It is easily verified that the coefficients satisfy the differential equations 
di = - igid 3 , d 2 = -ig 2 d 4 , 

<i 3 = 2zAirf 3 — igidi, d^ = 2iA 2 d 4 — ig 2 d 2 , (18) 

where 2Aj = (uq — ujj) is the detuning of the jth cavity frequency from the atomic 
transition frequency. 

Equations (1181) form decoupled pairs of differential equations that can be easily 
solved by using e.g. the Laplace transform technique. A simple solution of the equations, 
valid for an arbitrary initial state is given by 

di(t) = {fiidi(O) cos(ftit) - % [Aid^O) + gi d 3 (0)} sin^t)} , 

d 3 (t) = ^— {nxd 3 {0) cos(fiit) + i [Aid 3 (0) - ^idi(O)] sin(Q^)} , (19) 

where Vt\ = \/ g\ + A^ is a detuned Rabi frequency, di(0) and d 3 (0) are the initial 
values of the probability amplitudes. Corresponding expressions for the coefficients d 2 (t) 
and G?4(t), which determine the evolution of second cavity system, are obtained from (|T9l 
by exchanging 1 — >• 2 and 3 — > 4, and replacing gi by g 2 , and Ai by A 2 . 

The probability amplitudes oscillate sinusoidally with the Rabi frequency 
and their dynamics is strongly affected by the modulation term that depends on the 
detuning Ai and the coupling strength gj between the atom and the corresponding 
cavity mode. These modulation terms vanish periodically at time t n = nw/Qi. It is 
worth noting that the detuning enters the solutions in an antisymmetric way, whereas 
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the coupling strength enters the solutions in a symmetric way. This difference will be 
evident in the features of the time evolution of entanglement in the system. It is also 
worth pointing out that each of the amplitudes is influenced by the initial value of the 
other. This fact will also have a crucial effect on the evolution of an entanglement in 
the system. 

Solutions (fT9l) differ from previous results of Ref. [21] in that non-zero detunings 
and unequal coupling strengths are included. Non-vanishing detunings and differences 
between the coupling strengths lead to new aspects of the entanglement creation and 
evolution that differ qualitatively and quantitatively from those observed under exact 
resonance and equal couplings. 



3.2. Single- excitation entanglement concurrences 



Our interest is centered principally on the evolution of entanglement between the 
different parts of the two two-qubit subsystems 1 and 2. By denoting the two atoms as 
A and B, and the corresponding cavity modes as a and b, we may distinguish six pairs 
of sub-systems AB, ab, Aa, Ab, Ba and Bb, as investigated by Yonac et al [23l El]- 

To calculate the entanglement between any two qubit pair IJ, for example the 
atoms A and B, we take the trace over the other subsytems, to evaluate the reduced 
density matrix, from which the eigenvalues can be calculated. In general, the reduced 
states are two qubit mixtures, and have a density matrix of the general "X-form" [TTJIE] 



PU= n n n > ( 20 ) 
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which has the concurrence 

C u = 2max{0, \z\ - \fb~c}. (21) 

We solve for all the two qubit concurrences, Cab, C a b, Cao,, Csb, CUb, Csai and find 
that in terms of the probability amplitudes, their time evolutions are given by the simple 
expressions 

2|d 1 (t)||d 2 (t)|, C ab (t)=2\d 3 {t)\\d A (t)\ 



Cab® 

C Aa (t) = 2|d 1 (t)||d 3 (*)l, C aB (t) = 2\d 2 (t)\\d 3 (t)\ 
C Ab (t) = 2|d 1 (t)||d 4 (t)|, C Bb (t) = 2\d 2 (t)\\d 4 (t)\. 



(22) 



There are several remarks that can be made about the general results for the concurrence 
measures, before we proceed to the detailed analysis of the entanglement evolution. In 
the first place, we observe that an entanglement between the subsystems 1 and 2 cannot 
be created in time if initially the excitation was entirely in the 1 or 2 subsystem, so that 
Cab(£), C a f,(t), CAbit), Csaif) are all initially zero. Since the subsystems 1 and 2 are not 
directly coupled to each other, no transfer of the excitation, and hence entanglement, 
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is possible during the evolution. This means that an initial entanglement between the 
subsystems must be imposed to observe any entanglement evolution. 

Secondly, we note that an entanglement between A and a or between B and b can 
be created in time if the excitation is either in the 1 or 2 subsystem. This is easy to 
understand. Since the atoms interact with the cavity modes to which they are coupled, 
they can share the excitation with the corresponding mode. 

Finally, if initially C^a = 1 or (7b& = 1, to indicate maximum possible entanglement 
between each atom and its local cavity mode, one can easily show that the entanglements 
remain constant (frozen) in time unless the detunings Ai or A 2 are different from zero. 
In this case, the entanglement measures evolve in time, but cannot be zero. 

3.3. Initial states that result in a "frozen" pairwise entanglement evoultion 

There is the question of which initial states have stable, or "frozen", entanglement 
properties under the evolution of the J-C Hamiltonian. We consider the superposition 
with a uniform population distribution over the available energy states (EJ). We write 

l*o) = \ + e*|¥ a ) ± (|* 3 > - e^ 4 ») , (23) 

where 9 and <fi are arbitrary phase factors. For these states, in the case of zero detunings, 
we see from (1221) that all of the two-qubit entanglement between pairs remains constant 
with time. In this case, the entanglement is maximally shared (we call this "delocalized" ) 
between the six subsystem pairs, so Cab(O) = 0^(0) = .... = Cb&(0) = 1/2 (see Fig.|2J), 
in accordance with the rule (fl2|) . Thus, we see that where detunings are zero and 
for symmetric coupling constants gi — g2, the superposition state ( J23|) is a stationary 
state of the system, and hence must maintain fixed entanglement features. Interestingly 
though, this stability of all two-qubit entanglement concurrences is also predicted for 
asymmetric couplings, where g\ 7^ g 2 - 

The effect of detunings A on the time evolution of the atom-atom concurrence 
Cab(£) is quite different. For the initial state (T2"3"j) . the evolution of concurrence Cab 
is plotted in Fig. [2] for equal coupling constants g\ — g% — g. The atomic concurrence 
varies in time only for nonzero detuning and the time evolution of the concurrence is 
not symmetric with respect to the sign of the detunings. The dependence on time 
of the atom-atom concurrence CAB{t) changes qualitatively as well as quantitatively 
when the detunings A are varied from positive to negative values. A large and even 
maximal entanglement between the atoms can be created when the detunings are 
positive, whereas the initial atomic entanglement is reduced and can even be suppressed 
when the detunings are negative. We note that from the equality (fl2|) . a maximum 
entanglement (Cab = 1) between atoms will imply zero entanglement between the two 
cavities. When the atom- atom entanglement is maximum, at Cab — 1 ; the cavity- 
cavity entanglement is zero, at C a b = 0, and vice versa, so that this result also implies 
a channeling into the cavity-cavity system. In addition, an important result is that the 
maximum atomic entanglement Cab = 1 is created when the ratios Ai/pi = A 2 /g 2 = 1- 
Otherwise, the atomic entanglement is reduced. 



Entanglement evolution of two remote and non-identical Jaynes-Cummings atoms 9 



The effect of detuning implies that one can engineer the direction of evolution 
of entanglement by controlling the detuning. By this we mean that an initial "non- 
localized" fully symmetric entanglement and population, as given by ( J23|) can be 
transferred to a desired "localized" atom-atom entanglement and population by a 
suitable choice of the detunings: positive detuning can channel entanglement entirely 
into the atoms, or cavities, at appropiate times. "Localized" in this context means that 
the entanglement exists solely between the two atoms, or between the two fields, as 
in ©. 

3.4- Entanglement evolution for the Bell state \& 

Next, we consider the evolution of the concurrences for the following Bell initial state 

|tf> = (|*x) ± |* a » = (| U 00) ± I |T 00)) , (24) 

where the entanglement is restricted to just the two atoms. 

Here the population is solely localized in the atoms, and the entanglement of the 
pair is maximum, corresponding to Cab — 1- The time evolution of the concurrence 
measures is given by 

C AB {t) = |cos (Cligt) 1 1 cos (£l 2 gt) |, C ab (t) = |sin (Qigt) ||sin (Q 2 gt) |, 
C Aa (t) = |cos(fii5ft) ||sin(fii^) |, C Ab (t) = |cos (S^igrt) \\sm(tt 2 gt) \, 
C Ba (t) = \cos (Q 2 gt) \ \sm(Qigt) |, C Bb (t) = \sin (tt 2 gt) ||cos (£l 2 gt)\, (25) 

for the case of exact resonances A 1 = A 2 = and unequal coupling strengths, g\ ^ g 2 . 
Here fii = (1 + u/g) and Q 2 = (1 — u/g) are dimensionless Rabi frequencies, with 




Figure 2. Time evolution of the concurrence measure Cab showing stability of 
entanglement between atoms for zero detunings Ai = A2 = 0, and arbitrary g\, 
52, for the intial state ([25]) (solid line). The dashed lines reveal the effects of detuning, 
where g\ = 92 = g- Ai = A2 = g (dashed line), Ai = A2 = —g (dashed-dotted line), 
Ai = A 2 = 2.g (dotted line). 
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9 = (9i+ 92)12 and u={g l - g 2 )/2. 

For equal coupling strengths g\ = g<i = g and non-zero detunings, the bipartite 
concurrence measures are 



C AB {t) = cos 2 (Qfft) 
1 



sin 2 (Qgt) 



C ab (t) 



tt 2 



sin 2 (tit) 



C Aa {t) = C Ab (t) = C Ba {t) = C Bb {t) = - sin (Slgt) 



x 



cos 2 (Qgi) 



sin 2 (Clgt) 



1/2 



(26) 



where we have simplified the detunings to Ai = A2 = A and have introduced a 
dimensionless Rabi frequency, the same for both subsystems, Cl = (1 + 5 2 ) 1//2 with 
5 = A/g. 




Figure 3. Oscillation of the concurrence measures for equal coupling strengths. 
.9i = .92 = .9 and exact resonance Ai = A2 = 0, for the initial Bell state given 
by (©. C AB {t) (solid line), C ab {t) (dashed line), and C Aa [t), C aB (t), C Ab {t), C bB {t) 
(dashed-dotted line). Full transfer of entanglement from atoms to cavity modes takes 

place at t = n7r/2f2, where n = 1,3,5 , in accordance with both conservation 

rules (JTSJ) and |13]). Only ((T2j) holds for asymmetric g's (see Fig. (gj). 

These expressions are plotted in Fig. [3] to reveal the periodic transfer of 
entanglement from atoms to cavity modes, and vice versa. The oscillatory behaviour 
of the atom-atom concurrence was reported by Yonac et al. |24j . We note the transfer 
process of entanglement from the atoms to the cavity modes does not involve just the 
pairs C AB and C ab . By inspection of the time evolution of the concurrence measures in 
Fig. El we find that the initial maximal entanglement between the atoms is not only 
totally transferred to the cavity modes, but at the same time an additional pairwise 
entanglement is created during the evolution. This can also be seen by summing the pair 
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concurrence measures to find that at times t n = nir/AQ, where n = 1,3,5..., the total 
pair entanglement is larger than one. As time progresses, a part of the entanglement 
is transferred into the other pairs of the subsystems. Then, after a further interval, 
the entanglement is completely transferred into the cavity modes. The additional 
entanglement vanishes at times t n = nir/2Q, n = 1,2, 3, when the transfer process is 
completed. 



1 
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^ 0.6 
o 
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O 0.4 
0.2 



"0 1 2 3 4 5 6 
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Figure 4. Time evolution for the single-excitation state, showing the conservation 
of the sum of the individual pairwise 1 — 2 concurrences (SSPC): Cab (solid line); 
Cab (dottted line); C12 (dashed-dotted line); Cab + C a b (dashed line). The total 
c ab + C lb + C Ab + C Ib is ecmal t0 C i2 = 1 as indicated by l[T2|). Here A = 0, 
.9i = 2,92, 9 = (.91 + 52)/2 and a — 7r/4. We note that in this case, Cab + C a b is not 
constant, and the conservation rule (|13p breaks down. 

We note however that the nonlocal entanglement, as defined by (fT2|) is conserved 
throughout the transfer process. For the case of g\ = g 2 the rule ({TBI) also applies, as is 
evident in Fig. [3] Where there is asymmetry in the interaction, so we have non-identical 
cavities, this rule breaks down, though the general conservation rule (fl~2]) still applies. 
This is evident in Fig. [Hand has been noted by Chan et al. |27j . 

Another interesting prediction of (!25j) is that complete entanglement transfer from 
the atoms to the cavity fields requires exact resonances but not necessarily equal coupling 
constants. To examine the effect of detuning first, in Figures [5] and El we note that where 
detunings Ai and A 2 are nonzero, the localized initial entanglement is not completely 
transferred to another pair of qubits (see Fig. [7]). Detuning increases the oscillation 
frequency and decreases the minimum entanglement between the atoms. This is in 
complete contrast to the case of the initially completely symmetric entanglement, where 
entanglement transfer is possible only for a nonzero detuning. It is interesting, however, 
that the atomic entanglement returns to its initial maximum value periodically with 
the detuned Rabi frequency. This can be understood in terms of the localization of the 
initial energy. When A / 0, only a part of the initial energy, proportional to u is 
transferred to the cavity modes leaving the excess energy unlocalized. We illustrate this 
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situation in Fig. where we plot the concurrence measure Cab^) as a function of time 
and the detuning A = Ai = A 2 , and note that a nonzero detuning can thus be used to 
stabilise, or "freeze", the entanglement between the atoms. 

A second advantage to be given by imperfectly matched and asymmetric systems 
lies in the use of asymmetric couplings. We examine the effect of g\ 7^ <?2, which has in 
part been studied by Sainz and Bjork [TJ] and Cavalcanti et al [15]. These authors reveal 
that a periodic revival of full entanglement between atoms is possible, as well as a full 
transfer of entanglement to cavity systems. By studying all the pairwise concurrences, 
we report the new result, that for cavities with unequal rather then equal coupling 
strengths, one can engineer the transfer process so that the initial entanglement is fully 
transferred into a desired qubit pair. It is of practical importance to know if an initial 
localized entanglement can be transferred on demand with the perfect fidelity to a 
particular pair of qubits. Consider the time evolution of the initial atomic Bell state 
as given by (1241) . In the case of unequal coupling strengths, the entanglement evolution 
is described by (1251) from which one can easily find that if the ratio n g = g\/ g2 is not an 
integer number or a fraction of an integer number, no complete transfer of the state (1241) 
is possible to any of the qubit pairs. The complete transfer is possible only if n g is an 
integer number or a fraction of an integer number. However, the destination to where 
the initial entanglement can be completely transfered depends on whether n g is an even 
or an odd integer number. 

If n g is an even number, the initial maximal entanglement between the atoms, 
Cab(O) = 1, can be completely transferred only to the atom-field qubit pair C a b. On 
the other hand, if n g is a fraction of an even number, the initial entanglement between 
the atoms can be completely transferred only to the atomic- field qubit pair Csa- We 




Figure 5. Effect of detuning on the entanglement evolution of the Bell state "J, given 
by The concurrence measure Cab is plotted as a function of time and detuning A. 
Here g\ = g-2 = g. Detuning stabilises the entanglement between atoms and decreases 
the period of entanglement revival. 
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Figure 6. Effect of detuning on the evolution of the atom-atom concurrence Cab, 
for asymmetric <j's, and where the initial state is the Bell state ^ , given by (|24|) . Here 
gi = 2g2, and g = (gi + <?2)/2. Effect is similar to the case of symmetric <?'s, except 
for the modulation of the period of oscillation. 




Figure 7. Time evolution of the concurrences Cab (solid line), C a b (dotted line) for 
the Bell state showing the effect of detuning. The transfer of entanglement from 
atoms to cavities is not complete, but a full periodic transfer back to atoms occurs. Here, 
A = 2g and g\ = g? = g- We note that in this case, the simple conservation rule (fT5|) 
of Yonac et al [53] holds. 
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illustrate this situation in Fig. [HI where we plot the concurrence of the different qubit 
pairs as a function of time for exact resonances but unequal coupling strengths with 
n g = 2. It is evident from the figure that at particular discrete times, the initial 
entanglement between the atoms is completely transferred to the qubit pair Csa- The 
transfer of the entanglement to other pairs of qubits is incomplete at all times, the 
concurrences oscillate between zero and certain values below unity. The reason for this 
feature of the entanglement transfer can be understood intuitively by noting that for 
n g = 2 the Rabi frequency Oi of the population oscillation in the subsystem 1 is twice 
that of the Rabi frequency Jl 2 for the population oscillation in the subsystem 2. This 
means that over a complete Rabi cycle Q\t = it, the initial population in the subsystem 
1 returns to the atom, but at the same time the population makes a half Rabi cycle 
in the subsystem 2, i.e. the excitation in system 2 will be in the cavity mode. Thus, 
CAb = 1 at that time, with the concurrence in the other qubit pairs equal to zero. 

We now consider the situation where the ratio n g = g\j g2 is an odd integer number. 
Like the previous case of even n g , the initial entanglement can be completely transferred 
to only one of the qubit pairs, in this case, to the cavity modes. This is illustrated in 
Fig. [9j where we plot the time evolution of the concurrence of the qubit pairs for the 
exact resonance but unequal coupling strengths with n g = 3. During the evolution, 
the initial entanglement Cab{0) = 1 is completely transferred to the cavity modes, 
C a b{t) = 1 at Qit = 2tt. Again, the transfer of the entanglement to only one of the 
entanglement pairs can be understood in terms of the Rabi cycles of the population in 
each of the subsystems. Since Qi is triple Q2, over the one and half of a Rabi- cycle 
in the subsystem 1, the population is completely transferred to the cavity mode a, and 




Figure 8. Effect of asymmetric coupling constants on the evolution of the 
entanglement of the initial Bell state 4", given by Concurrence measures are 

evaluated for resonance Ai = A2 = and gi = 2g2'- Cab (i) (solid line), CAb(t) (slim 
line). In this case, maximum entanglement is transferred periodically from atoms AB 
to the atom-field system Ab. 
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Figure 9. Effect of asymmetric coupling constants on the evolution of the 
entanglement of the initial Bell state 'J', given by (|24|) . The concurrence measures 
are evaluated for Ai = A2 = and g\ = 3g2- Cab (t) (solid line), C a b(t) (dashed line). 
In this case, maximum entanglement is transferred to the field-field system ab. 



at the same time, the population in the subsystem 2 is also completely transferred to 
the cavity mode b. Consequently, at this time, the cavity modes become maximally 
entangled with C a b(t) = 1. 

In concluding this section, we would like to emphasize that we have found a 
significant difference between the evolutions of initial "fully symmetric non-localized" 
and "targeted localized" entanglements. That is, we have found that the "non-localized" 
entanglement of states (J23l) evolves only if the atoms are imperfectly coupled to the 
cavity modes, so that detunings A are nonzero. On the other hand, where entanglement 
is initially targeted solely to the atoms, as in fT2"4"|) . large detunings provide a means of 
stabilising the entanglement so it cannot be completely transferred to the field system. 
Moreover, for zero detunings in this case, one can achieve complete transfer of an initial 
entanglement to a desired pair of qubits, this transfer process being steered by changing 
the ratio between the coupling strengths g of the atoms to the cavity modes. 

4. Evolution of two-excitation entanglement 

We now turn to the problem of entanglement dynamics when two excitations are can 
be present. This allows us to consider the evolution for the case of the Bell state 

|$) = cosa| tt 00) + e ie sina\ || 00), (27) 

which has been studied previously [23], and found to undergo an abrupt depletion of 
entanglement for a 7^ 7r/4, called "sudden death of entanglement" (ESD). 

Using the Schrdinger equation, we set up equations of motion for the probability 
amplitudes and then solve them to study the dynamics of the entanglement. It is easy 
to show that the equations of motion for the probability amplitudes of the states ((8]) 
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group into three independent sets of coupled equations. We define amplitudes for the 
wavef unction as 

|#> = d 1 {t)\ rr oo) + d 2 (t)\ n 01) + d 3 (t)\ it io) + d 4 (t)\ II ii) 

+ do(*)|UOO) (28) 

in terms of the energy basis states defined by (jSJ). The evolution of this set involves four 
equations 

d 1 = - i(Ai + A 2 )d~i - i(g 2 d 2 + 01^3), 
d 2 = - i(Ai - A 2 )d 2 - i(g 2 d~i + gid*), 
d 3 = i(Ai - A 2 )d 3 - i(gidi + g 2 d<i), 

d 4 = i(A t + A 2 )d 4 - i( gi d 2 + g 2 d 3 ). (29) 

We note the remaining equations that would arise by including the possibility of two 
quanta per site form two independent sets, namely 

= —i\Plg\d§, d~e = iAide — i\Plg\d^, (30) 

and 

d 7 = -iV2g 2 d 8 , d 8 = iA 2 d 8 - i\^2g 2 d 7 , (31) 

where states 5 — 8 are | || 10), | II 20), | It 01), | II 02) respectively. We have removed 
the fast oscillating terms by introducing a rotating frame through the relation 

= e'^+^di, i = l,2,...,8. (32) 

As before for the single excitation, we have allowed for the possibility of imperfect 
couplings of the atoms to the cavity modes by introducing detunings Ax and A 2 , and 
unequal coupling strengths g\ and g 2 . 

Note that the dynamics of the coupled amplitudes d 5 (t) and d 6 (t) as well as d 7 (t) 
and ds(t) involve, in fact, an exchange of a single photon between the atoms and the 
cavity modes to which they are coupled. They can display an entanglement during the 
evolution, but the dynamics will be similar to that discussed in Sec. Ill for a single 
excitation. Therefore, we will not consider the evolution of these states and will focus 
instead only on the evolution of the other two-photon states that involve simultaneous 
evolution of two excitations. 

Although the set of equations of motion (|29|) involves four coupled equations, it is 
analytically solvable. We can put Eqs. (!29l into a matrix form and solve them by the 
matrix inversion. Due to the complexity of the general solution, we study separately 
the time evolution of the probability amplitudes for two cases: (1) non-zero detunings 
and equal coupling strengths, (2) zero detunings and unequal coupling strengths. 

Instead of the four equations for di (i = 1,2,3,4), we introduce the following 
combinations of the probability amplitudes 

d«(t) = di{t) + d A {t), d {2 \t) = - d A {t), 

d {3) (t) = d 2 (t) + d 3 (t) , d (4) (t) = d 2 (t)-d 3 {t), (33) 
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and find that in the first case, the time evolution of the amplitudes is given by 
d^(t) = rf {1) (0) cos(Qgt) - ^ [d (3) (0) + <M (2) (0)] sin(^t), 

d?>(t) = ^ [S 2 \0) - 5d^(0)} + ^ [S 3 \0) + 8dP(0)] oos(n^) 

-— ^(O)sin(^t), 

S 3 \t) = ^ [5S 3 \o) - S 2 \o)} + A [d* 3 >(o) + ^(o)] oos(n^) 

-|^(0)sin(^t), 

rf( 4 )(t) =rf( 4 )(0), (34) 

where the detunings are simplified to Ai = A2 = A, so we have introduced a 
dimensionless Rabi frequency, the same for both subsystems, f2 = 2(1 + 5 2 ) 1 / 2 with S = 

For the second case of exact resonance but unequal coupling strengths, we find 

d^(t) = d {1) {0) cos{2gt) - id {3) (0) sm(2gt), 
d®(t) = d {2) (0) cos(2wt) + id iA) {0) wn(2ut), 
d( 3 \t) = d {3) {0) cos(2gt) - id {1) (0) wn(2fft), 

d^ (t) = d {4) (0) cos(2wt) + id {2) (0) sin(2wt) , (35) 
where, as before for one-photon states, g = (gi + g<i)j2 and u = (gi — g<i)j2. 

4-1. Entanglement evolution for the double- excitation state 

We first investigate the evolution of the bipartite entanglement with two quanta present 
in the system, but where there is always one excitation in each of the subsystems 1 
and 2. The initial state of the system is a linear superposition of the form 

|tf> = di(t)|TT 00) + d 2 (t)\U 01) + d 3 (t)\ |T 10) + d 4 (t)\ ii 11). (36) 

This state cannot give any entanglement in the atom-atom subsytem, regardless of the 
values of the coefficients involved. This is readily seen by evaluating the reduced density 
operator for the atom systems. We evaluate 

Pab = Tr cavity p = ^ (i\(j\p\j)\i), (37) 
m=o,i 

where \z) and \j) refer to cavity modes a and b respectively, to get 

pab = Mi(*)| 2 | TT) + W)\ 2 \ Tl) + \d 3 (t)\ 2 \ t|) + \d,(t)\ 2 \ |l), (38) 

which is a mixture of separable states and hence cannot be entangled. The atom-atom 
bipartite entanglement is always zero, CAB{t) = 0. The same conclusion applies to 
the bipartite entanglements C a f,(t), C^t) and Cs a (t). The only non-zero concurrence 
possible is for the qubit pairs CAaif) and Csbif)- In other words, an entanglement can 
be created during the evolution between the atoms and the cavity modes to which they 
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are coupled. The amount of entanglement that can be created in these pairs can be 
determined from the relations 

C Aa (t) = 2\d 1 (t)d ¥ 3 (t) + d 2 (t)dl(t)\, (39) 

and 

C Bb (t) = 2|dx(t)d*(t) + d 3 (t)dl(t)\. (40) 

Figure [H]] shows the entanglement creation and evolution between the atom A and the 
cavity mode a, with an initial separable state | f 00). Similar entanglement properties 
appear between the atom B and the cavity mode h. We see that an entanglement 
between the atom and the cavity mode is created during the transfer process of the 
excitation between them. The maximum entanglement is observed at times when the 
excitation is equally shared between the atom and the cavity mode, and vanishes when 
the excitation is completely located at either the atom or the cavity mode. In terms of 
the Rabi oscillations, the system becomes disentangled at every half of the Rabi cycle 
of the oscillation, i.e. at gt = nn/Q (n = 0, 1, 2, . . .). 




Figure 10. Evolution of the concurrence CA a (t) for the initial separable state | "ft 00) 
and two different detunings: A = (solid line) and A = g (dashed line). 

Similar to the case of single excitation states of Section [3j the system becomes 
entangled for a longer time when the frequencies of the atom and the cavity field are 
detuned from each other. In other words, an imperfect matching between the atom 
and the cavity field leads to a more stable entanglement than in the case of the perfect 
matching. In terms of the Rabi oscillations, the system becomes disentangled only at 
every Rabi cycle of the oscillation, i.e. at gt = 2nii/Q. This is easy to understand 
if one considers entanglement as resulting from a superposition corresponding to a 
derealization of energy. When the frequencies of the atom and the cavity mode are 
different, the initial amount of energy localized in the atom is not completely transferred 
to the cavity mode. A part of the energy, not absorbed by the cavity mode remains 
delocalized, which results in a nonzero entanglement. 
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4-2. Inclusion of the auxiliary state |0) 

The above analysis has showed that no pairwise entanglement is possible between any 
parts of the two subsystems 1 and 2 when we have with certainty a single excitation in 
each of 1 and 2. We can introduce entanglement by including in the wave function an 
auxiliary state, the ground state |0) for which there is no excitation [23J: 

|0) = |U00). (41) 

The evolution of the entire wave function can be written explicitly as 

|* > = di(t)|TT oo) + d 2 (t)\ u oi) + d 3 (t)UT io) + <k(t)\u n) 

+ do(*)|U 00). (42) 

It is easy to show that the state ( T4T1) is not the eigenstate of the Hamiltonian (CQ). 
Therefore, the probability amplitude d does not evolve in time, d (t) = d (0). As before 
for the single-excitation states, we will examine the entanglement between different 
pairs of the subsystems and its transient behaviour especially under imperfect matching 
conditions. For this general case, the pairwise concurrences between the subsystems are 
found to be 

C AB (t) = max{0, 2(|d 1 (*)||d (*)l - \d 2 (t)\\d 3 (t)\)} , 
C ab (t) = max{0,2(\d 4 (t)\\d (t)\-\d 2 (t)\\d 3 (t)\)}, 
C Ab (t) = max{0,2(|d 2 (t)||doW|-MiWI|d4WI)}, 

C aB (t) = max{0,2(|d 3 (t)||do(*)| " \di(t)\\d,(t)\)} ■ (43) 

Evidently, we see that entanglement between the pairs of subsystems is possible only if 
the ground state is included. 

We note interesting properties of the concurrence and differences between the 
single and double-excitation cases. Firstly, the concurrences appear as differences of 
products of the absolute values of the probability amplitudes. This gives a possibility 
for discontinuities in the behavior of entanglement that was absent in the case of the 
one-excitation, where the concurrences, given by (122j) . are continuous functions. The 
abrupt disappearance of entanglement that persists for a period of time is referred to 
as "sudden death of entanglement" (ESD) and was reported by Yonac et al. [23] for 
this case. Secondly, there is evident a competition in the creation of entanglement 
Cab between two pairs of states (|$i), \&4)) and (($2), l^))- Entanglement creation in 
concurrence pairs involving the states and |<&4) is diminished by the presence of 
population in the states |<J>2) and and vice versa creation of entanglement involving 
the states |$ 2 ) and |$ 3 ) is diminished by the presence of population in the states |$i) 
and I $4). In terms of the population transfer, entanglement creation by a simultaneous 
exchange of two photons is diminished by one-photon exchange processes, and vice versa. 
The competition between these one and two-photon processes has been recognized as 
the source of the phenomenon of ESD, sudden death of entanglement. 
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4-3. Inequality for sum of square of the nonlocal pairwise concurrences: 

The reduced density matrix for the atom-atom system, in terms of basis states | ft), 
IT!), |!T) and HI), is 
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dependence for simplicity of notation. 




Figure 11. Evolution of entanglement for the Bell state <& showing the lack of 
conservation of the sum of the square of the nonlocal pairwise 1 — 2 concurrences, 
SSPC = \Cab\ 2 + \C Ab \ 2 + \C aB \ 2 + \C ab \ 2 . Here A = and gi = g 2 = g. The squares 
of concurrences are plotted in bold for a = n/6, and as thin lines, for a = w/12. 
We have C 2 AB (solid line), C 2 ah (dashed line), C\ h = C 2 B (dash-dotted line). As the 
initial entanglement weakens, to a = n/12 for example, regions of total loss of all non- 
local pairwise entanglement are evident (SSPC=0). Total pairwise 1 — 2 entanglement 
(SSPC) is restored when the entanglement has been completely transferred from atoms 
to cavity modes. 

In contrast with the single-excitation case of Section [3J there is the difference that 
the reduced state can be entirely separable, for periods of time (sudden death), even 
when the off-diagonal coherence term d\d^ which is the source of the entanglement is 
non-zero. This feature results in a reduction of the sum total of the squares of the nonlocal 
pairwise concurrences (SSPC) and is also responsible for the entanglement sudden death 
(ESD) phenomenum. The entanglement for the remaining reduced systems, the two 
cavities, and atom and cavity at different location, is given similarly, by the concurrences 
as in (1431) . 

In fact, for the Bell state $ of the form ([6]), we can identify regions of evolution 
for which each of Cab, C a b, C a B, C a b is zero, meaning that all nonlocal pairwise 
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entanglement involving systems 1 and 2 is lost. This effect has been discussed by 
Chan et al [27], and is evident in Fig. [HI for a = n/12. In this case, the entanglement 
involves the four states, that of the atoms and cavity modes. We readily see from ( 1431) 
that where we have \di\ = \d 2 \ = \d-s\ = = |, for any nonzero value of do = sin a, 
each pairwise concurrence between systems 1 and 2 will be zero. The entangled state 
at t = nir/AQ, where n = 1,3,5... , that can give this total depletion of all pairwise 
entanglement for tana < 1/2, is written 

|cosa(| TT 00) + HI 01) + lit 10) + 11)) + sina|!! 00). (45) 
4-4- Controlling ESD using detunings 

Sudden death of entanglement has been reported by Yonac et al. [23J, for the case where 
the initial state is the Bell state $ 

|$) = cosa| TT 00) + sina| II 00), (46) 

and where the system has zero detunings and equal coupling constants. They showed 
that for a ^ 7r/4, entanglement is lost in a discontinuous fashion, as displayed in Fig. 
[T2l for a = 7r/12. As tana deviates from 1, the length of time for which entanglement 
is lost is increased. 

We observe from the full solutions that an imperfect matching of the cavities 
to the atomic transition can prevent the onset of ESD, to provide a stabilization of 
entanglement. This effect, which is similar to that of Figure [U is shown in Fig. [T2l 
and may be understood by realizing that the nonzero detuning A inhibits the full 
entanglement transfer from atom to cavity modes. 




Figure 12. Plot of the evolution of Cab for the intial Bell state <!>, given by ((6|). Here 
.9 = .9i = 92 and a = 7r/12. The effect of detuning A ^ is to remove the ESD, the 
sudden death of entanglement. 
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For some initial states however, we point out that the detuning will act to induce 
ESD, that is, to destabilise the entanglement, though stabilization is reached for 
sufficiently large A. To further investigate the role of the imperfect matching on the 
evolution of entanglement, we consider two different initial states, one that displays ESD 
with zero detuning A, and another that displays stable entanglement with A = 0. We 
consider a linear superposition of the zero-excitation state with a two-quanta symmetric 
state 



i n oo) + m ii» + m oo) 



(47) 



and a linear superposition of the zero excitation state with a two-quanta antisymmetric 
state 
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Figure 13. Concurrence between the atoms as a function of time t and detuning A for 
two different initial states: (a) the asymmetric state (|35|) and (b) symmetric state (jTT)) . 



Figure [13] shows the evolution of concurrence Cab(£) as a function of the 
detuning A = Ai = A 2 of the cavity modes from the atomic resonances. We see 
from the figure that the detuning has an opposite effect on the two initial states: the 
stable entanglement at A = for the antisymmetric state is forced into an entanglement 
sudden death (ESD) with the onset of a non-zero A, whereas the unstable entanglement 
evolution (ESD) displayed by the symmetric state at A = is stabilized by detunings 
A / 0. We note in both cases however, that for very large A, the entanglement is 
stabilized at or near the intial value. 
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4-5. Effect of asymmetric coupling constants 

The loss of pairwise entanglement for the spin correlated Bell state $ is a feature 
maintained even in the case of asymmetric couplings: g\ 7^ gi- The depletion 
of entanglement (ESD) for the atom-atom pair has been reported, for asymmetric 
couplings, by Sainz and Bjork [TJ]. In Fig. [TH we plot the non-local pairwise 
concurrences from which one can calculate the sum of their squares, SSPC, as in (JUJ), 
to reveal a loss of total pairwise entanglement where couplings are asymmetric. As with 
the symmetric case, the effect is more pronounced with a decreasing initial entanglement, 
Ci2- We note from Fig. [15] that there is a shift of the period of entanglement sudden 
death, as the ratio of the g's increases. 




Figure 14. Evolution of the nonlocal pairwise concurrences for the Bell state $, where 
couplings are asymmetric: gi = 2gi, A = 0. Here a = ir/6. The entanglement sudden 
death (ESD) is evident, as is the lack of conservation of the sum of the square of the 
pairwise 1-2 concurrences, SSPC = \Cab\ 2 + |CUb| 2 + |C a B| 2 + |C b| 2 . Individual 
concurrences are: C\ B (solid line), C 2 h (dashed line), C% B (dash-dotted line), C\ h 
(dotted line). Also plotted is C\ 2 > the horizontal dotted line at 0.75. 

In Fig. [TB], we plot the atom-atom concurrence Cab{^)i for different coupling 
strengths, for the different initial state 

|A> = -^(| IT 10} + 111 00)), (49) 

which has no entanglement between the two atoms. In contrast to the state $, there is 
no abrupt depletion (ESD) of entanglement. We note that a suitable choice of coupling 
ratio can lead to a periodic maximal entanglement bewteen the atoms. 

In summary of this section, we point out that small differences in the frequencies can 
lead to the disappearance of the zeros (ESD) in the entanglement evolution, whereas 
small differences in the coupling strengths do not generally prevent the vanishing of 
entanglement, but rather tend to create ESD. In other words, entanglement sudden 
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death is much more sensitive to an imperfection in the frequencies than to the coupling 
strengths. 

5. Conclusion 

We have studied the time evolution of entanglement between two remote and in general 
non-identical Jaynes-Cummings systems. A particular attention has been paid to the 
effect of an imperfect matching of the atoms to the cavity fields and to the understanding 
of the nature of entanglement evolution and steered transfer. We have quantified 
entanglement through the pairwise concurrence between different parts of the two 
systems and have found that entanglement transfer can be triggered and controlled by a 
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Figure 15. Evolution of atom-atom concurrence for the Bell state $, with a = tt/6, 
where couplings are asymmetric: A = and g\ — <?2 (solid line); <?2 = 2gi (dashed 
line); g2 = lOgi (dashed-dotted line). The ESD shifts toward a later time as the ratio 
of coupling constants increases. 




Figure 16. Evolution of the concurrence CUs(i) for the initial state |A) with g\ = gi 
(solid line), gi = 2g\ (dashed line) and gi = 0.5gi (dashed-dotted line). 
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suitable choice of the detunings between the atomic and the cavity field frequencies and 
the coupling strengths of the atoms to the cavity fields. We consider two cases: the first 
where there is only one excitation present in the entire system; the second where there 
is a superposition of two excitation levels, zero excitation and two excitations, one in 
each of the two remote systems. These two cases exhibit completely different properties 
of the entanglement evolution and steering of the entanglement transfer. When the 
atomic and field frequencies are detuned from each other and there is only one excitation 
present, an initial entanglement cannot be completely transferred between the pairwise 
subsystems. On the other hand, for zero detunings one can achieve complete transfer of 
the entanglement to a desired pair of the subsystems and this process can be steered by a 
suitable choice of the ratio between the coupling strengths. In the second case where two 
excitations are present in the system, we have showed that no entanglement is possible 
between any parts of the two subsystems when the number of photons in each subsystem 
is conserved. An additional auxiliary state of zero excitation has to be included to create 
entanglement between the subsystems. We also discuss a conservation rule governing the 
pairwise entanglement between the systems, that the sum of the squares of the pairwise 
concurrences is conserved during the evolution of an initial entanglement, for the one- 
excitation case. For the superposition of the zero and two excitation case, the sum is 
not conserved, and can be zero for intervals of time, implying an entanglement sudden 
death in all pairwise concurrences. Finally, we have addressed the issue of sudden death 
of entanglement in relation to non-identical cavities, and have found that a mismatch 
in the frequencies can lead to the disappearance of sudden death in the entanglement 
evolution, whereas differences in the coupling strengths can prevent and even create the 
vanishing of entanglement. 
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